Quasiparticle anisotropy and pseudogap formation from the weak-coupling 

renormalization group point of view 
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Using the one-loop functional renormalization group technique we evaluate the self-energy in the 
weak-coupling regime of the 2D t-t' Hubbard model. At van Hove (vH) band fillings and at low 
temperatures the quasiparticle weight along the Fermi surface (FS) continuously vanishes on ap- 
proaching the (n, 0) point where the quasiparticle concept is invalid. Away from vH band fillings 
the quasiparticle peak is formed inside an anisotropic pseudogap and the self-energy has the con- 
ventional Fermi-liquid characteristics near the Fermi level. The spectral weight of the quasiparticle 
features is reduced on parts of the FS between the near vicinity of hot spots and the FS points 
closest to (n,0) and (0,7r). 
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The apparent non-Fermi liquid (nFL) normal-state fea- 
tures and the intriguing pseudogap structures remain 
a key issue in the ongoing research activities on high- 
T c superconductors. Angle resolved photoemission spec- 
troscopy (ARPES) keeps providing new and increasingly 
detailed insight into the nature of the correlated elec- 
tronic states responsible for the unconventional normal 
state properties and the superconducting pairing state 
[1,2]. Refined ARPES data continue to serve as chal- 
lenging test grounds for theoretical concepts developed 
for studying correlated lattice electrons in two dimen- 
sions (2D). 

The pseudogap phenomenon in cuprates manifests it- 
self by the absence of well-defined quasiparticles (qps) 
near the (tt, 0) point of the Brillouin zone (BZ) and the 
seeming partial disappearance of the FS [2]. One di- 
rection of recent research has been to seek the origin 
of these apparent non-FL characteristics and pseudogap 
structures in the competing nature of different ordering 
tendencies indicative for the existence of quantum critical 
points in the ground-state phase diagram [3,4]. 

On the other hand, the vicinity to the transition from 
a metal to an antiferromagnetic (AF) insulator in doped 
cuprates demands the understanding of the role of AF 
spin fluctuations on the evolution of spectral properties 
on the metallic side of the metal-insulator transition. 
Anisotropic spectral features naturally arise within the 
2D Hubbard model near its AF instability even in the 
weak-coupling regime [5]; they were studied within the 
spin-fcrmion model [6], fluctuation exchange [7,8] and 
two-particle self-consistent (TPSC) approximations [9] 
and more recently cluster extensions of the dynamical 
mean-field theory (DMFT) [10] and cluster perturbation 
theory [11]. QMC studies on finite clusters [12] provided 
valuable insight into the doping evolution of the spec- 
tral function on a sparse momentum grid in the strong 
to intermediate coupling regime, but they were not able 
to trace the pseudogap formation into the weak-coupling 



regime. The results of previous studies agree with each 
other only for the case of perfect nesting (next-nearest 
neighbor hopping t' = and half-filling) where a pseudo- 
gap opens everywhere around the FS upon cooling. The 
experimentally most relevant case of a non-nested FS re- 
mains an issue of debate. 

From the apparent necessity for the development of al- 
ternative methods for the analysis of correlated electron 
systems in 2D, functional renormalization group (fRG) 
techniques have proven to be most promising additional 
tools with an emphasis on the identification of the leading 
electronic instabilities in the weak-coupling regime of the 
2D Hubbard model and its extensions [13-17]. Although 
the electronic on-shell lifetime, the FS shift and the ap- 
proximate qp weight calculated from the imaginary part 
of the self-energy at the lowest Matsubara frequency were 
recently considered within the many-patch fRG analysis 
[18], these quantities alone are only indicative but not 
sufficient to trace the evolution of the spectral properties 
in the vicinity of the AF instability. 

In this letter we use the one-loop fRG scheme for the 
t-t' weak-coupling Hubbard model to evaluate the self- 
energy on the real frequency axis. We show that at vH 
band fillings the qp weight along the FS continuously van- 
ishes from the BZ diagonal towards the (tt, 0) point where 
the qp concept is invalid. The qp weight suppression is 
accompanied by the growth of two additional incoherent 
peaks in the spectral function, from which an anisotropic 
pseudogap originates. On moving away from vH band 
fillings, a qp peak with small weight is formed inside the 
pseudogap. This peak regains spectral weight from the 
incoherent pseudogap features on moving towards the BZ 
diagonal. These results offer a complementary weak cou- 
pling view of pseudogap formation from the RG perspec- 
tive with new implications for the doped AF-insulator 
regime of the Hubbard model. 

Specifically, we consider the Hubbard model for N e 
electrons on a square lattice 
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FIG. 1. Imaginary part of the self-energy at imaginary 
frequencies within the one- loop RG approach at U = 2t, 
t' /t = 0.1 and vH band filling n=0.92. ki? is chosen in the 
patch closest to the vH point at (ir,0). The dots show the 
positions of the first and second Matsubara frequency. 
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where the hopping amplitude tjj = t for nearest neigh- 
bor (nn) sites i and j and tij = —t' for next-nn sites 
(i, t' > 0); for convenience we have shifted the chemical 
potential /i by At' . We apply the fRG approach for one- 
particle-irreducible functions with a sharp momentum 
cutoff [15]. This technique considers the effective action 
obtained by integrating out modes with energy \e-k\ > A 
where = — It (cos k x + cos k y ) + At' (cos k x cos k y + 1 ) — 
is the electronic dispersion; A (0 < A < A = max|ek|) 
is the cutoff parameter. This procedure is used in the 
weak-coupling regime for small and intermediate t' < 0.3i 
when the ferromagnetic instability is absent [16]. In this 
scheme the flow of the self-energy Ea(1c, iuS) is given by 



dA 



= V a °Sa 



(2) 



where o is a short notation for the summation over 
momentum-, frequency- and spin-variables, see e.g. Ref. 
[15]. The renormalization of the electron-electron inter- 
action vertex V\ at one-loop order is given by 



dV A 
dA 



V\ o (G a o Sa + Sa o Ga) o Va- 
The propagators Ga and 5a are defined by 
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(4) 



We neglect the influence of the self-energy and the FS 
shift on the RG flow and therefore the self-energy is not 
included in the Green functions (4). Eqs. (2) and (3) 
have to be solved with the initial conditions Va = U 
and Ea = 0. Since also the frequency dependence of the 
vertices is neglected, it is convenient to reinsert, following 
Ref. [18], the vertex from Eq. (3) into Eq. (2) to obtain 
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FIG. 2. Self-energy in SOPT (dashed line) and one-loop 
RG (solid line) at U = 2t, t'/t = 0.1, and vH band filling 
n=0.92 at T = 0.17t (Vmax = Wt). kj, = (3.100,0.003) is 
chosen in the first patch, closest to the vHs. The insets show 
the self-energy and the spectral function at small frequencies. 
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To solve Eqs. (3) and (5) numerically we divide the mo- 
mentum space into 48 patches with the same patching 
scheme as in Re fs. [15,16]. The self-energy on the real axis 
is obtained by analytical continuation using Pade approx- 
imants [19]. To resolve fine structures close to the Fermi 
level, we take advantage of Eq. (5) that for frequency- 
independent vertices the self-energy can be calculated at 
arbitrary frequencies on the imaginary axis. A sample 
result for kp in the first FS patch closest to the (tt, 0) 
point is shown in Fig. 1. Obviously, the imaginary part 
of the self-energy behaves non-monotonically at frequen- 
cies a; < u>i = irT with [d Im£(kir, itu) / 801}^^+ > 0. By 
the Cauchy-Riemann conditions for an analytic function 
this necessarily implies [<9Re£(kF, oj) / dbj] u ^ M + > and 
thus the breakdown of the qp concept at this particular 
value of ki?. Evaluating the self-energy at the Matsubara 
frequencies iu) n only (circles in Fig. 1) this behavior is 
missed. Therefore, for the analytical continuation we use 
a dense set of frequencies on the imaginary axis at small 
|cj| < t and a set of iui n for \u>\ 3> t. This procedure can be 
viewed as an analytical continuation first from Matsub- 
ara frequencies ioj n to some suitably chosen set of points 
on the imaginary axis and a subsequent continuation to 
the real axis using Pade approximants. The quality of 
the approximants was checked by both the analysis of 
the analyticity in the upper half-plane [19] and by re- 
quiring the fulfillment of the sum rule for the resulting 
Green function. 

First we consider the results at the vH band filling 
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FIG. 3. Same as Fig. 2 for T = 0.082* (V max = WOt). 
In addition, c) and d) show the spectral functions and the 
qp weights for kj- in different patches around the FS. The 
spectral functions are shifted by (n — l)0.5t for a better view; 
for n = 1, 2 they are multiplied by 6. 

(jz = 0) for t' = O.lt and U = 2t. The self-energy at tem- 
perature T = 0.17f calculated within fRG together with 
the result in second-order perturbation theory (SOPT, 
which is obtained by the replacement V — > U in Eq. 
(5)) is shown in Fig. 2. Remarkably, in the first patch 
closest to the vH singularity (vHs) at (tt, 0), SOPT and 
fRG results show both a minimum of ImE(kir,cj) at the 
Fermi level u> = instead of a maximum as expected for a 
FL. Simultaneously, ReS(ki?, to) has a positive slope near 
u> = 0. The dip in ImS calculated within fRG is much 
more pronounced than in SOPT. Moreover, the physical 
origin of these features is very different. The peculiarities 
of the SOPT self-energy arise solely from the thermal ex- 
citation of electrons near the vHs and although they exist 
at all T > 0, they lose their relevance with decreasing T. 
Instead, the more pronounced anomalies in the fRG self- 
energy increase in size with decreasing T and reflect the 
growing AF correlations and the concomitant tendency 
towards pseudogap formation. Importantly, the single 
peak in A(kp ,u) — — ImG(k^,Ci; + i0 + )/-7r in Fig. 2 is 
not a qp feature because the low energy structure of £ 
invalidates the qp concept. 

In Fig. 3 we decrease the temperature to T = 0.082i 
which is closer to the crossover temperature T* where the 
one-loop interaction vertices V tend to diverge [15] as a 
fingerprint of the nearby AF instability. Although the 
low-temperature regime where the effective vertices are 
so large (V max = lOOi) is outside the validity region of the 
one-loop approach, it proves most useful to clearly iden- 
tify the structure of the spectral function which originates 
from the above discussed form of the self-energy. In the 
first patch a two-peak structure in A(kp ,u>) arises with 



a local minimum at the Fermi energy. Outside the first 
patch, the real part of the self-energy has a narrow region 
with negative slope near uo — accompanied by a local 
maximum of ImS (Fig. 3b), which leads to the formation 
of a qp peak inside a pseudogap structure (Fig. 3c). This 
behavior is in fact reminiscent of the finite temperature 
DMFT results in the metal to insulator crossover region 
[20]. For larger \u>\ the behavior is qualitatively similar 
in all patches. The qp peak quickly merges with the two 
incoherent peaks of the pseudogap on approaching the 
BZ diagonal along the FS. The qp weight around the FS 
(see Fig. 3d) gradually vanishes with approaching the 
(tt, 0) point. For t' = and the corresponding vH band 
filling n = 1, the qp peak is absent in all patches. All 
spectral functions along the FS show in this case a two- 
peak pseudogap structure at low temperatures similar to 
the TPSC results [9] and the dynamical cluster approach 
[10]. 

In Figs. 4 a,b we show the self-energy for U = 2t, 
t' = O.lt at high and low temperatures, respectively, for 
/i = 0.03t (n — 0.94) when the Fermi level is slightly 
above the vHs energy. At higher temperatures only weak 
anomalies in S(k^,w) are present near the Fermi level. 
Upon lowering the temperature towards T* a clear maxi- 
mum of ImS develops at the Fermi energy similar to that 
in Im£(kj>,a;) for n > 1 at vH band fillings. The real 
part of self-energy has the proper negative slope for a FL 
in a narrow energy window \uj\ < /i around the Fermi 
level. In the spectral function (see Fig. 4c) we observe 
the split-off of two incoherent peaks near the qp peak at 
the Fermi level. With increasing \x and the corresponding 
increase of temperature which is kept above T* , the qp 
peak gains weight from the incoherent pseudogap peaks. 
For [i < (fillings below vH band filling) the value of 
T* rapidly drops with \/j,\ for finite t', and we observe a 
quick crossover to a single-peak structure. 

Fig. 4d shows the anisotropic weight of the qp feature 
in the 3-peak spectral functions of Fig. 4c for different 
chemical potentials /i > 0. Away from the vH band filling 
the qp weight is finite in the first patch near the (tt, 0) 
point, which is now below the Fermi level. Hot spots, 
i.e. points on the FS which are connected by Q = (tt, tt) 
and thus most affected by the scattering from AF spin 
fluctuations, are identified as local shallow minima in the 
angular dependence of the qp Z-factor along the FS line. 
The central qp peak of the spectral function is suppressed 
on parts of the FS between the near vicinity of hot spots 
and the FS points closest to (tt, 0) and (tt, 0) and quickly 
gains spectral weight upon moving towards the BZ diag- 
onal. 

In summary, we have investigated the self-energy on 
the real frequency axis in the 2D t-t' Hubbard model 
within a one-loop fRG analysis. For vH band fillings the 
self-energy has a non-FL form at the FS points (tt, 0) and 
(0, tt) which are connected by the AF wavevector Q. Qps 
exist everywhere else on the FS but with anisotropic spec- 
tral weight. The RG flow indicates that at low tempera- 
tures the continuous decrease of the qp weight along the 
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FIG. 4. (a,b) Self-energy for a patch closest to the vHs at 
U = 2t, t'/t = 0.1, /i = 0.03i, T = 0.18t (V^oa = Wt; a) 
and T = 0.093t (V moai = 100*; b). (c,d) The evolution of the 
spectral function and the quasiparticle weight for different 
chemical potentials, the temperatures are chosen such that 
Vmax=W0t. The lines in (c) are shifted by 0.5i for a better 
view. 

FS is accompanied by the simultaneous growth of two ad- 
ditional incoherent peaks in the spectral function. Away 
from vH fillings a qp peak at (n, 0) and (0, it) emerges 
with small spectral weight. 

These results provide a complimentary weak-coupling 
scenario for pseudogap formation and offer implications 
for the evolution of anisotropic spectral properties at fi- 
nite temperatures on the metallic side of the transition 
to an AF insulator. The strong angular dependence of 
the qp properties in our fRG study is reminiscent of 
the recent high-resolution ARPES studies on underdoped 
cuprates [2]. Although the fRG method does not allow 
to go below the crossover temperature T* , and therefore 
to explore the spectral functions at a fixed temperature 
in a wide range of dopings, we expect further suppres- 
sion of the weight of the central peak with decreasing 
temperature and/or increasing interaction strength, and 
simultaneously the increase of the doping range near half 
filling where the pseudogap structures emerge. 

The main change of the low-energy physics upon mov- 
ing from the weak- to the strong-coupling regime is ex- 
pected in the redistribution of spectral weight from the 
vicinity of the Fermi energy to the developing lower- and 
upper Hubbard subbands. These subbands should form 
in addition to the low-energy structures we have eluci- 
dated in our fRG analysis. This conjecture finds indeed 
support from recent results of the dynamical cluster ap- 
proximation [10] and cluster perturbation theory [11]. 



The observed features are distinctly different from the 
previously found low-temperature crossover to a non-FL 
form of the self-energy at half filling [7,9] as well as from 
the proposed hole pocket picture and the partial destruc- 
tion of the FS in strong coupling Hubbard or t-J models 
[21]. The fRG results thereby offer a valid alternative for 
the interpretation of experimental data from photoemis- 
sion spectroscopy on cuprates. 
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